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Abstract. An interaction of a photon with 3 H is invstigated based on a three dimensional Faddeev approach. 
In this approach the three-nucleon Faddeev equations with two-nucleon interactions are formulated with consid- 
eration of the magnitude of the vector Jacobi momenta and the angle between them with the inclusion of the 
spin-isospin quantum numbers, without employing a partial wave decomposition. In this formulation the two 
body t-matrices and triton wave function are calculated in the three dimensional approach using AVI 8 potential. 
In the first step we use the standard single nucleon current in this article. 



1 Introduction 

Since the early days of the study of the nuclear physics so 
many efforts have been performed on 3N systems consider- 
ing real or virtual photon interactions [1 1-[2|. Also several 
studies on the behavior of 3N bound states in real or virtual 
photon absorbtion have been reported[3|-[4]. Before six- 
ties variational approach was used for these calculations 
and works using this approach are still continuing. After 
introducing Faddeev formulation for three body systems 
0-161, new efforts using this scheme were started. As an 
example one can point out the early calculations of elec- 
trodisintegration Q and photodisintegration|8| with 3 He 
and 3 H. An improvement in the photodisintegration calcu- 
lation of the bound and 3N continuum with the same 3N 
hamiltonian have been performed [9 |. There are also other 
approaches to calculate electromagnetic interactions with 
light nuclei such as Green-function-Monte-Carlo method 
IfTOl . hyperspherical harmonic expansions IfTTIl . and Lorentz 
integral transform (LIT) method [ 12]. There is a very good 
review of Faddeev calculations on the interaction of real 
or virtual photon with 3 He [13|. In this work like previ- 
ous calculations the partial wave decomposition has been 
used. In PW approach one should sum all PW to max- 
imum angular momentum where the calculation is con- 
verged. The problem is that in higher energies this max- 
imum angular momentum increases and we should solve 
more complicated equations. To avoid this complexity one 
should use vector momentum as basis states [14|. To this 
aim in the past decade the main steps have been taken 
by Ohio-Bochum collaboration (Elster, Glockle et al.) and 
Bayegan et al. to implement the 3D approach in few-body 
bound and scattering calculations (see for examples Refs. 
irT5l -122l). It should be clear that the building blocks to the 
few-body calculations without angular momentum decom- 
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position are two-body off-shell t-matrices, which depend 
on the magnitudes of the initial and final Jacobi momenta 
and the angle between them. Fachruddin et al. have cal- 
culated the NN bound and scattering states in a 3D repre- 
sentation using the Bonn-B and the AVI 8 potentials lfl"5ll - 
fl6l . Recently there has been efforts to do the same calcu- 
lation using chiral potential [23 1. our aim in this work is 
to formulate photodisintegration of 3 H in a three dimen- 
sional Faddeev approach. In the first step we ignore three 
body forces and we just use the single nucleon current. We 
will use AVI 8 potential and triton wave function which has 
been calculated in our previous work EDI . 

This manuscript has been organized as follow: in sec- 
tion |2] we explain our basic states and we evaluate all of 
the matrix elements in these basis. In section [3] we intro- 
duce our singularity problem and its solution. We finish in 
section |4] with a summary and outlook. 



2 Integral equation of nuclear matrix 
elements without partial wave 
decomposition 

To calculate the photodisintegration observable we first need 
to calculate nuclear matrix elements in the Faddeev scheme. 
For more details see Ref. lfT3l . 



N = |<0ol(l + tG )P\U) 



\U) = (1 + P)J\t/f) + tG P\U) 



(1) 



(2) 



In above equations t is NN t-operator which obeys Lipmann- 
Schwinger equations, Go is free propagator, P is permuta- 
tion operator, is three body bound state and \U) is an 
axillary state. Three body forces have been ignored. |0o) is 
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a subsection of the fully antisymmetric free state, |<?o)> m 
which nucleons 2 and 3 are in subsystem. 



|3>o> = (l+P)k*o> 



(3) 



|0o) is a ls° our basic state to solve the integral equation 
(2) and is antisymmetric under permutation of nucleons 2 
and 3. 



|0o> = |pqmim 2 m 3 viV2V3) a 



(4) 



In equation (4) p and q are jacobi momenta and m's 
and v's are the spin and isospin of the individual nucleons 
respectively. 

Orthonormality and completeness relations of these ba- 
sic states can be considered as bellow: 



a {pqm 1 m2m3ViV2Vi\p'q'm' l m' 2 m' 3 v' l v' 2 v' 3 ) a 

-6(p + p')£m 2 m' v'jWq-qO'Vm'Aiv; (5) 



q|pqmim 2 m 3 viV2V3) a 



m\ m 2 m 3 

V\ V 2 V 3 



x fl <pqm 1 m 2 m 3 v 1 V2V3l = 1 (6) 



Considering these properties we can rewrite the integral 
equations (1) and (2) in our basic stats. 

N = ^ a (wmm 2 m3ViV2V3\(l +tG )P\U) 



1 



"<pqmim 2 m 3 viv 2 y 3 |P|f/> 



+ - a {pqmim 2 m2ViV2V3\tG P\U) 



(7) 



fl (pqmim 2 m 3 ViV 2 v 3 |t/> 

= a (pqmim2m3Viv 2 v3|(l + P)J\\fi) 

+ a (pqm l m 2 m 3 v l v 2 v 3 \tGoP\U) (8) 

The effect of permutation operator on our basic states 
can be considered as follow: 

(pqmim 2 m3ViV 2 v 3 |P|p'q'm' 1 m 2 OT 3 v 'i v 2 v 3 ) 
= <*(P + \v' + ^q'Wq - p' + ^q') 

^&m\m 2 Om 2 m' 3 ^m^m' x $v\ v' 2 v' 3 v' { 

1 3 1 
+S(p + -p' - -q')6(q + p' + -q') 

= 6(p + n 2 )6(p' - tti) 



Where: 

^i=q+^q' n2=\q + q' (10) 

Now with respect to above relation and symmetry consid- 
erations we can evaluate equations (7) and(8) as follow: 

113 1 

N = 2^~2 P ~ 4 q ' P ~ 2 qm2mmV2V3Vl ^ 

1 3 1 
(-^P + 4*1' ~P ~ -qm3mim2V3V\V2\U)} 



^ J d 3 q a (p m 2 m 3 v 2 v 3 \t\^q + q', m 2 m 3 v 2 v 3 > fl 



m 2 m 3 



q2 + 1 q ,2 +q .q, <~^Q' ~ Q. q' m'^m^V^U) (11) 



<pq, m 1 m 2 m 3 v 1 V2V3|f/> 

= <pq, m 1 m 2 m3Viv 2 V3l(l + P)J\0 

+ ^ J"^V<P OT 2«3V2V3|f|^q + q^ 2 m 3 v 2 v 3 ) fl 
q77q^7q-^7<-^q'-q>q'«»iv 2 v 3 vi|f/) (12) 

^ m 

The firs term in the equation(12) can be evaluated as: 
fl <pq, mim 2 m 3 v 1 v 2 v 3 |(l + P)J\\fi) 
= J] f rfVrfV 

m',v' ^ 

fl (pq, mim 2 m 3 viV2V 3 |(l + P)7|p'q', m^m^y^)" 
x fl <p'q', m'jm^m^v'jv^v^) (13) 

Now we concentrate on the elements of these equations i.e. 
current, two-body t-matrix and triton wave function, more 
precisely. 

2.1 current 

Considering the symmetry properties we have: 

fl <pq, miOT 2 m 3 viv 2 y 3 |(l + P)J SN \ifs) 

= 3 fl <pq, mim 2 m 3 vi v 2 v 3 |(l + P)J SN (l)\i//) 



(14) 

Matrix elements of single nucleon current can be evaluated 
as follow: 

fl (pq, m 1 m 2 m 3 viv 2 v 3 |7(l)|p'q', m\m' 2 m' 3 v\v' 2 v' 3 ) 
= <5(q'-q+^Q) 

X 2^(P ~~ P')dm 2 m 2 dm 3 m' 3 d V2 v 2 Sy iV > 3 

~ <^(P + P')^m 2 m^m 3 m^V2V^V3V^] X ^ (Q, q) 



Vl V 



(9) 



(15) 



19 International IUPAP Conference on Few-Body Problems in Physics 



In above equation Q is the momentum of photon. We need 
to rewrite the single nucleon current operator in a form 
which is suitable for our basic states. The current opera- 
tor which we will use is: 

J = G E (Q)^-^ + -^-G M (Q)crx(k l -k' l ) (16) 

2mpj zmjv 

Which is summation of convection current and spin cur- 
rent. Ge(Q) and Gm(Q) are electric and magnetic form fac- 
tors respectively. For the convection part we have: 

k, +k; = 2q + Q+^K (17) 

As we will show we have to choose coordinate system in 
which the z axis is along the Q vector and we also need 
tensor component of current so the second and the third 
terms of the right hand side of equation(17) will vanish. 
Thus for the convection current we have: 

jcon.ec = Ge(Q) ^ (18) 

m N 

And the tensor component of spin part can also be eval- 
uated as: 



f A ^(p',p, C ose) = v^(p',p,e) 
i 

+ \ J dp"p" 2 J d(cose")v^ A (p',p",0\6")G o (p") 
-l 

t? A '(p",P,0") 

1 

+ \ J dp"p" 2 J d(cos8")v* A %' A (p\p",e',6")G (p") 
-l 

C( P ",P,8") (22) 
Where 

In 

vf,' A A (p\ P ",e',6") = J ^V^-^C ( P ',p") 



(23) 

2.3 triton wave function 



IniM 



For evaluating the Triton wave function we need to make 

(19) 

a relation between this wave function in our basic states to 
the one which has been calculated in the following basis 1 20 



2.2 two-body t- matrix 

Two body t-matrices can be related to the one which cal- 
culated in helicity basis: 

fl <pm 1 OT 2 nv I V2|f|p'miOT2v' 1 V2> a = ^(v 1 +v 2 ),(y I +i/ 2 ) 

71 St 

C( - -5 , mim 2 A )C( - i S, m' l m' 2 A' ) 

y JN((t, p -<p' )jS id \jS iqi \ 



AA 
nSl 



d A'A^°PP'^ 



f A %(p,p' 7 COS0 pp ,,z) 



(20) 



In the above relation z — E - 22- is the energy of subsys- 
tem. As we know two-body function has a singularity in 
the energy of deuteron, z = Ed- To remove this singularity 
we should consider t-operator as follow: 



a (pmim2nv\V2\t\p'm\m2v' 1 v 2 ) a 
a (pmim2«Vi V2\i\p' mim2v' l v' 2 ) a 
Z-E d 



(21) 



Two body t-matrix in helicity basis has been calculated 
beforeHa. 



(M(^)Sm s (t-)Tm T m = (pq,X pq ,a\^) (24) 



We can relate these states to our free spin and isospin states 
with Clebsch-Gordan coefficients. 



Where 



g 7 a = (y\a) 



\a) = \(s l -)Sm s ,{t\)Tm T ) 
\y) = \mim 2 m 3 viV2V 3 ) 



(25) 



(26) 



It is very important to mention that the spin of the nucle- 
ons is quantized in the direction of the z axis which in the 
calculation of wave function it has been chosen to be in the 
direction of q. But we have to consider the z axis along the 
direction of incident photon Q. So we should first rotate 
the spin of the nucleons in our basis to be settled in the 
direction of q axis. Then we should use Clebsch-Gordan 
coefficients to obtain the wave function in the calculated 
basis mentioned in the equation (24): 

<pqm i m 2 /«3 vi V2V3 llA} 

= ^ ^^D mm ' i (8 q ,ip q )D mim >JJ) q ,(p q )D mm 'J,9 q ,(pq)gy a 



x(pq,X pq ,a\tlr) 



(27) 
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3 Singularity problem 

In order to consider the singularity problem we can rewrit- 
ten the equation (11) an (12) in a unified form ignoring 
isospin dependent which is similar to spin dependent. 



(p,q,Q) = f/; 

I illzlll} (p,q,Q) 

U, 



Zf ,3 ,, m 2 m 3 



fc,q",Q)C»«;-.; ( P' jr| ' z) 



(28) 



To solve this integral equation we should evaluate singular- 
ity in the denominator of the propagator which is a func- 
tion of q" and angle between g"and q. So instead of sin- 
gular point we have a region of singularity in q — q" plane. 
There is a solution to this moving singularity in Ref.[24|. 
For using this method we have to put z axis along the q. 
But because of simplification in current operator and final 
cross section we should choose the z axis in the direction 
of the momentum of the photon, Q. So in order to evaluate 
the singularity we should use another method which is in- 
troduced in Ref.|25|. Therefore one should separate angle 
part of delta functions as follow: 

5(p +7ii)6(p -ni)= — 

<5(p'+^)5(p"-^i) 
6( P ' - m) 6(p" - ax) 

p' 1 p" 2 
tf^'-^W+rfi) (29) 



<5(p'-7Ti)(5(p"+;ri) 



And then the integral equation can be rewrite as follow: 

(p,q,Q) = t4 (p,q,Q) 

+ ^ f d 3 q"dp'dp" 

U m ' m ' mi (7T2, q"> Q) ^m 2 m 3 m'm'. ;(P> 711 ' Z ^ 



E-h(P" 2 + h" 2 )E + ie-E d - 3 ^ 



4m 



(30) 



After some simplification the integral equation transforms 
to this equation: 

11112111} (p,q,Q) 

CO 

2 



1 



q/2+p' 



mi mi ~ m VJ 4 1 ' 

J" dq"q"G(q,q",p') j dq"6(x" - x Q ) 

\ql2-p'\ 

U m ^ mi (p"7T 2 , q", Q)C 2B!3m> ; (P< P% , z) 

CO 

m 2 m 3 



9/2+9" 

^ dp' p'G(q,q" , p') J dq"6(x" - xq) 

19/2-9" I 

l/^m, (p"*2, q", Q)C«^ (P« ^> • z > (3 1 > 

In the above equation G which is always positive is defined 
as: 



G(q,q",p') = 



1 



-E d - 3 -£ + h(p' 2 + h 2 ) 



(32) 



x" = cos 0" indicates the angle between q and q" and 
xq is introduced as follow: 



1 / /2 1 2 //2\ 1 / "2 1 «2 2\ 

x =—(p --.q-q ) = —{p ~ ~.q -q) 
qq" 4 gg" 4 



(33) 



4 Summary and outlook 

In this paper we have formulated the Faddeev integral equa- 
tions for calculating the photodisintegration observable of 
triton in a three dimensional approach. To this aim we in- 
troduced our basic states which contains jacobi momenta 
in vector forms as well as individual spin and isospin of 
each nucleon. So we have avoided to decompose angle 
states in terms of angular momentum states (partial wave 
approach) which is traditionally used to solve these kind of 
equations. The final integral equations are less complicated 
than the PW ones and are unique in number of the equa- 
tions in all energies. We have also explained about over- 
coming of the moving singularity in our work. 

The calculation of this observable using the AVI 8 po- 
tential is underway and the results will be published soon. 

Adding two and three body currents as well as three 
body forces in our calculations are other future major works. 
The same calculation for radiative capture is also under 
consideration. 
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